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SUMMARY 


Our  work  during  the  past  year  on  active  ELF/VLF  generation 
has  produced  major  breakthroughs  in  the  understanding  of 
nonlinear  interactions  in  a  dissipating  medium  such  as  the 
ionosphere  for  the  purpose  of  downconverting  H.F.  to  ELF/VLF 
waves.  A  key  result  was  the  fact  that  contrary  to  the 
conventional  understanding  of  such  processes,  the  presence  of 
a  dissipative  medium  lowers  the  threshold  for  complete 
decay  and  increases  the  downconversion  efficiency  over  the 
one  expected  on  the  basis  of  the  Manley  Rowe  relations. 

Two  codes  were  developed  and  optimized. 

(i)  The  first  code  named  CMPOWER  requires  as  input  a  model 
of  the  ionospheric  electron  density  and  magnetic  field, 
and  the  transmitter  characteristics.  It's  output  is 
the  ELF/VLF  signal  as  a  function  of  the  modulation 
frequency  and  distance  from  the  HF  facility. 

(ii)  The  second  named  TWO  (three  wave  coupling)  is  a 
nonlinear  code,  which  includes  dissipation  and  finite 
pump  size  effects. 


I.  INTRODUCTION 


Generating  VLP/ELP  waves  by  utilizing  the  ionosphere  as 
an  active  medium  is  an  exciting  prospect,  since  it  avoids 
many  of  the  difficulties  connected  with  large  and  inefficient 
ground  based  VLP  and  ELP  transmitters.  By  adjusting  the  HP 
beam  geometry,  the  low  frequency  waves  can  be  made  to  propa¬ 
gate  upwards  in  the  magnetosphere  for  use  in  active  magneto- 
spheric  stimulation  experiments,  or  downwards  into  the  earth 
ionosphere  waveguide  for  communication  purposes.  Two  basic 
downconversion  mechanisms  have  been  proposed: 

(i)  Modulation  of  the  auroral  or  electrojet  current  by 
large  HF  heaters1*6. 

(ii)  Stimulated  nonlinear  coupling  of  two  HF  wavepackets 
with  frequency  difference  equal  to  the  desired  VLF/ 
ELF  frequency  4'7'8. 


The  first  mechanism  has  a  long  history,  dating  from  the 
"Luxembourg"  effect.  The  mechanism  is  based  on  the  possibil¬ 
ity  of  changing  the  ionospheric  electron  temperature  by  HF 
radiowave  heating.  The  electron  neutral  collision  frequency 
ven  dePends  almost  linearly  on  the  electron  temperature  Te, 
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where  the  "o"  subscript  denotes  properties  of  the  ambient 
ionosphere.  Consequently,  the  electron  conductivity  of  the 
lower  ionosphere  can  be  made  to  vary  in  time  at  a  frequency  w 


as  a  result  of  the  temperature  changes  caused  by  the  ampli¬ 
tude  variation  of  the  HP  waves.  The  "Luxembourg"  effect 
arises  because  the  time  dependent  conductivity  of  the  lower 
ionosphere  causes  a  time  dependent  attenuation  of  HF  radio- 
waves  propagating  through  it.  At  the  receiving  station,  this 
appears  as  an  amplitude  modulation  of  the  HP  waves.  The 
VLP/ELF  generation  by  current  modulation  is  based  on  the 
formation  of  time  dependent  current  loops  in  the  lower  iono¬ 
sphere  by  modulating  the  electron  temperature  in  regions  of 
naturally  occurring  geophysical  currents,  such  as  the  auroral 
and  the  equatorial  electrojets.  Note  that  current  loops, 
rather  than  electric  dipoles  are  expected  at  VLF/ELF  fre¬ 
quencies  due  to  the  large  conductivity  of  the  ionosphere 
(i.e.,  a  >>  oj/4h).  This  approach  has  been  very  successfully 
applied  in  the  Soviet  Union3  and  more  recently  to  the  MPE 
experiments  at  Tromso,  Norway1'2  and  by  the  Pennsylvania 
State  University  group  at  Aracibo9.  More  details  concerning 
this  approach  can  be  found  in  the  review  given  by  Stubbe10. 

The  second  mechanism  is  more  novel7'8  and  generates  the 
nonlinear  current  sources  without  relying  on  pre-existing 
ionospheric  currents.  It  mixes  two  HF  wavepackets  in  the 
ionosphere  to  produce  VLF/ELF  signals.  However,  unless  the 
process  occurs  in  a  stimulated  (i.e.,  lasing)  fashion,  the 
downconversion  efficiencies  are  trivial.  Stimulated  down- 
conversion  occurs  for  pump  power  above  a  threshold  value  that 
depends  on  the  ionospheric  parameters,  the  HF  and  VLF/ELF 
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frequencies  and  mainly  on  the  angle  of  HF  beam  to  the  iono¬ 
spheric  density  gradient.  Recent  theoretical  calculations 
indicate  fairly  high  downconversion  efficiencies.  An  impor¬ 
tant  variation  of  the  scheme  using  only  one  HF  pump  can 
amplify  pre-existing  ionospheric  VLF/ELF  wavepackets  by  tens 
of  db's.  This  mechanism  has  not  yet  been  experimentally 
proven,  although  evidence  of  the  amplifier  operation  has 
fortuitously  been  seen  in  Aracibo19. 

Our  work  in  the  past  year  emphasized  the  second.  In  the 
section  we  review  some  of  the  fundamental  physical  notions 
involved  in  collisionless  interactions.  Section  III  dis¬ 
cusses  some  recent  results  concerning  the  semi-collisional 
regime,  where  the  frequency  of  the  VLF/ELF  wave  is  below  the 
electron  neutral  collision  frequency.  This  regime  is  of 
interest  for  ELF  generation  in  the  E-region.  The  final 
section  applies  our  theoretical  results  to  the  ionospheric 
VLF/ELF  generation.  Appendix  I  discusses  the  basic  deriva¬ 
tion  of  the  equations,  while  Appendices  II,  III  include 
recent  publications. 

II.  MAXIMUM  ENERGY  TRANSFER  IN  COLLISIONLESS  THREE  WAVE 

INTERACTIONS 

The  equations  describing  the  generation  of  a  lower  fre¬ 
quency  wave  (a>3,k3),  from  the  interaction  of  two  high  fre¬ 
quency  wave  packets  (w1,k1)  and  (<*>2,k2)  satisfying  the 
resonance  conditions 

w  1  ~  w2  *  <*>  3 

(3) 

*1  “  *2  *  3 
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where  Uj  are  the  group  velocities,  Tj  the  damping  rates, 
and  aj  the  slowly  varying  complex  amplitudes  of 
packets,  normalized  to  action  density13'14,  so  that 


wave- 


a?  = 


E  • 


p- 
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(5) 


Ej  is  the  electric  field  amplitude  of  the  wave  and  Dj  is 
its  dielectric  function.  In  deriving  Equations  (4),  we  have 
assumed  positive  energy  waves  in  a  collisionless  ( oj1  >u)2>u)3>> 
rl'r2»r3)»  homogeneous  medium,  with  the  finite  size  of  the 
pump  (<*>1,^)  being  the  only  source  of  inhomogeneity  (Figure 
1).  The  value  of  <  for  any  conservative  interaction  of 
interest  can  be  found  by  following  the  method  of  electro¬ 
dynamics  of  weakly  non-linear  media  or  the  ponderomotive 
force  approach 1 1 ' 1 3  * 1 8 .  An  example  of  such  a  calculation  can 
be  found  in  Ref.  18. 

The  important  issue  in  the  downconversion  problem  is  the 
maximum  efficiency  of  power  transfer  from  the  pump  to  the 
lower  frequency  wave  and  the  conditions  under  which  it  can  be 
achieved.  For  the  collisionless  case  described  here,  the 
maximum  power  transfer  from  the  pump  Pj  to  the  low  frequency 
wave  P3  can  be  found  from  the  conservation  of  the  system's 


action  and  energy12  and  is  given  by  the  well  known  Manley- 
Rowe  relations  as 


(6) 


An  important  practical  issue  is  how  close  to  the  maximum 
theoretical  efficiency  given  by  (6)  can  one  reach  and  under 
what  conditions.  An  extremely  important  result  given  by 
Reiman18  by  solving  numerically  the  system  of  Equations  (4), 
is  that  if  we  define  a  length 


IU2,U3«l1/; 

| <ax ( t=0 ) j  ' 


(7) 


maximum  efficiency  (i.e.,  complete  pump  decay)  can  be 
achieved  if 

U2xU3x  <  0  (Sal 


where  L  is  the  size  of  the  pump.  Equations  (8)  a'e  strictly 
valid  for  a  square  pump,  but  they  have  been  found  to  vary  by 
less  than  a  factor  of  two  for  other  pump  profiles.  Figure  2 
shows  the  efficiency  n  achieved  as  a  function  of  L/Lc 
compared  to  the  theoretical  efficiency  given  by  Equation 
(6)  as  found  in  Reiman14,  demonstrating  that  h/ho  a  0.7  - 
0 . 8 n0 .  Note  that  the  numerical  results  are  consistent  with 
analytical  results  based  on  the  method  of  inverse  scat¬ 
tering14'16,  which,  however,  is  applicable  only  for  non¬ 
overlapping  initial  wavepackets. 


An  observation  with  profound  consequences,  is  that  the 
conditions  for  complete  decay  given  by  Equations  (8)  coincide 
with  the  conditions  for  an  absolute  instability 1 7  of  the 
system  of  linear  equations 

(IE  +  U2X  lx)a2  =  VX)a3 

(9) 

(|t  +  u3x  Ix^3  =  YQ(x)a2  . 


Equations  (9)  can  be  found  from  Equations  (4)  by  assuminq  the 
pump  energy  (i.e.,  ja^2)  as  constant. 

We  found  above  the  conditions  for  pump  depletion  for  the 
case  of  one  dimensional  space-time  evolution.  Of  more 
interest  for  our  application  is  the  steady  state  three  dimen¬ 
sional  boundary  value  problem.  If  the  group  velocities  of 
the  waves  are  coplanar,  the  steady  state  equations  reduce  to 


(U1z  ll  +  U1x  Ix^l  =  <a2a3 


^U2z  Iz  +  U2x  lx  5*2  ■ 


(U3z  h  +  U3x  l^a3 
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* 

<a  ^  a  2 


(  10) 


It  is  easy  to  see  that  the  results  of  the  previous  section 
are  directly  transferable  to  Equations  (10)  if  we  make  the 
transformation 
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Equations  (10)  then  become 
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Note  that  the  role  of  the  time  variable  is  now  played  by  z, 
so  that  the  pump  depletes  while  it  propagates.  Therefore  the 
steady  state  interaction  corresponds  to  an  amplifier  in  space 
for  a3,  while  the  absolute  instability  found  before  would 
correspond  to  an  oscillator  in  time.  The  conditions  for 
complete  decay  in  the  steady  state  interaction  would  be 


U2U3  <  0  (13a) 


L 


> 


ji 

1 


(13b) 


where 


a 

L 

c 


lu2u3j 1/2 

j  <  a3  (z=0) I 


(  14) 


For  U.  ,  Uj  ,  U,  >  0,  as  is  our  case,  (13a)  coincides  with 
(8a)  and  is  automatically  satisfied.  Using  the  definitions 
given  by  Equations  (11),  it  is  easy  to  see  that  the  threshold 
condition  (13b)  is  identical  to  (8b),  but  with  the  boundary 
condition  ai(z=0)  replacing  the  initial  condition  a^t^O). 
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In  summary,  the  results  presented  above  indicate  that 
downconversion  efficiences  of  the  order  of  70  to  80%  of  the 
theoretical  maximum,  i.e., 

P3  <*>3 

T|  «  -5-  -  (.7-. 8)  —  (  15) 

P1  “>1 

can  be  achieved  for  interactions  whose  daugther  waves  satisfy 
(8a)  and  pump  power  and  size  satisfy  (8b).  Note  that  the 
threshold  computations  involve  only  the  group  velocities  of 
the  daughter  waves  and  the  characteristics  of  the  pump  (i.e., 
P3  and  antenna  gain). 

III.  MAXIMUM  ENERGY  TRANSFER  IN  THE  SEMI-COLLISIONAL  REGIME 

The  analysis  presented  above  is  valid  in  the  collision¬ 
less  regime.  In  examining  the  generation  of  ELF  waves  in  the 
lower  ionosphere  by  HF  pumps,  we  notice  that  the  low  fre¬ 
quency  u>3  can  be  much  smaller  than  the  electron  neutral 
collision  frequency  ven  which  controls  the  damping  of  the 
HF  waves.  The  damping  of  ELF  waves  is  determined  by  the  ion 
neutral  collision  where  w3  >>  v^n.  We  call  the  regime 

“>i,  w2  >>  ven  >>  w3  >>  vin 

the  semi-collisional  regime.  In  this  regime  the  dominant 
part  of  the  low  frequency  nonlinear  force  is  due  to  the 
thermal  nonlinearity  rather  than  the  usual  ponderomot ive 
force.  This  leads  to  nonsymmetric  coupling  coefficients  in 
Equations  (4)  with  rather  profound  consequences.  The  coupled 
equations  for  the  semi-collisional  regime,  whose  derivation 
is  given  in  Ref.  18  are 


(■ST  +  -1  *  ■3*  +  r!  )ai 

+  -2  *  lx  +  r2^a2 

(^t  +  -3  *  lx  +  r3 ^a3 


(  .  .  4  v  \  * 

*( 1+17  z; Jaia2 


(  16a) 

(  16b) 

(  16c) 


where  <  is  the  coupling  coefficient  of  the  collisionless 
interaction  and  v  =  ven.  For  v  +  0,  we  recover  Equations 
(4). 


The  conditions  for  absolute  instability  can  be  found  as 
previously  by  assuming  the  pump  a!  constant  and  solving  the 
eigenvalue  problem  given  by  Equations  (16b)  and  (16c)  with 
the  appropriate  boundary  conditions.  We  find  that  Equations 
(8)  are  now  modified  to 
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For  v  <<  u>3,  (17b)  would  reduce  to  (8b).  For  v  >>  u>3, 
however,  the  pump  power  threshold  is  reduced  by  a  factor  of 
the  order  v/w3.  Since  the  system  of  Equations  (16)  is  not 
conservative,  we  cannot  apply  the  inverse  scattering  method 
to  the  three  wave  system.  To  determine  the  downconversion 
efficiency  when  the  conditions  (17)  are  satisfied,  we  solved 
Equations  (16)  numerically.  For  simplicity,  we  considered  a 
model  problem  of  a  one  dimensional  space-time  evolution  with 
U 1  x=  0,  U2x=  ~U3x,  r j  *  0,  =  10.1,  U)2  *  10.0,  V  =  1.2  and 
w3  =  .1;  jaj  =  1,  |a2|  =  .01,  |  a  3  j  =  0..  The  results  are 
shown  in  Figure  3.  The  interaction  leads  to  pump  depletion. 
More  importantly,  the  downconversion  efficiently  is  increased 
by  a  factor  of  the  order  v/w3,  over  the  value  given  by  (6). 
The  preferential  transfer  of  action  to  a3  is  understandable 
since  for  v  >>  w3,  the  coupling  coefficient  of  the  low  fre¬ 
quency  equation  is  v/ui3  larger  than  that  of  the  high  fre¬ 
quency  sideband.  Physically  it  means  that  additional  w3 
photons  are  generated  at  a  large  rate  via  a  temperature 
perturbation,  with  a  net  loss  of  energy  to  the  medium.  As  a 
comparison,  we  show  also  the  results  of  the  collisionless 
interaction.  In  order  to  have  the  proper  threshold 
normalization  for  a  .given  pump  width  L,  we  have  used  | a  x |  = 
1.0  for  the  collisionless  case  and  ja^  =  |A|~1/2  =  .25  for 
the  semi-collisional  case  so  that  the  respective  normalized 
widths  L/Lc  and  L/Ac  have  the  same  value.  From  Fiqure  3, 
we  see  that  pump  depletion  is  seen  to  occur  in  both  inter¬ 
actions.  For  the  collisionless  case,  the  emergence  of 
symmetric  pulses  identifiable  as  solitons  is  as  predicted  by 
the  1ST  solutions,  and  ja3j  *  |a2j  which  implies  downconver¬ 
sion  efficiency  of  the  order  of  u>3/w1  consistent  with  the 
Manley-Rowe  relations.  In  the  collisional  case,  however, 
|a3j  *  ( v/u3 ) 1 1  2 j a2 |  indicating  a  major  improvement  in  the 
efficiency. 
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IV.  APPLICATION  TO  IONOSPHERIC  VLP/ELF  GENERATION 

In  this  section  we  discuss  the  HP  antenna  requirements  in 
order  to  achieve  efficient  downconversion  from  HF  to  VLP/ELF 
in  the  ionosphere.  Let  us  first  summarize  the  highliqhts  of 
the  previous  sections. 

(i)  Our  study  was  confined  to  interactions  produced  by 
two  HF  pumps  propagating  transverse  to  the 
ionospheric  density  gradient.  Such  interactions 
minimize  the  required  threshold  for  complete  pump 
decay  by  avoiding  the  phase  mismatch  caused  by  the 
density  gradient.  The  practical  implications  of 
this  to  system  design  are: 

(a)  Beam  steerability  such  that  grazing  HF  inci¬ 
dence  can  be  achieved  at  ionospheric  altitudes 
of  optimum  interaction. 

(b)  Antenna  gain  such  that  the  beam  width  L  in  the 
direction  of  the  density  gradient,  is  at  most 
equal  to  the  density  gradient  scalelenqth  L^ 
( i.e. ,  L<Lj^) . 

(c)  Ability  to  radiate  several  frequencies  simu¬ 
ltaneously. 

(ii)  For  collisionless  interactions  (w3>ven)'  com¬ 
plete  pump  decay  can  be  achieved  for 


L  % 

rc  * 
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in  which  case  downconversion  efficiency  is 


-  (.7-. 8)  ~  . 


(iii)  For  semi-col  1  is  ional  interactions  ( w1 ,  u2>  ven>  w3 ) » 
complete  pump  decay  can  be  achieved  for 


>  1.8 


in  which  case  the  downconversion  efficiency  is 

p3  w3  ven 

■=—  -  (.7-. 8) - —  . 

P  ^  u  ^  oo  3 

We  have  applied  the  collisionless  interaction  to  VLF 
generation  at  an  altitude  of  150  Km  and  the  semi-col 1 isional 
interaction  of  ELF  generation  at  a  lower  altitude,  -100  Km, 
where  electron  neutral  collisions  are  important.  In  both 
cases,  we  considered  a  5  MHz  HF  pump  incident  perpendicular 
to  the  ionospheric  density  gradient  and  almost  parallel  to 
the  magnetic  field  (angle  between  the  wavevector  k^  and  B0  is 
5s)  with  a  beamwidth  L  equal  to  the  local  density  scalelength 
Ln«  We  have  taken  Ljq  to  be  20  Km  at  150  Km  and  10  Km  at 
100  Km.  The  threshold  electric  field  required  to  achieve 
pump  decay  was  computed  from  Equations  (8b)  and  (17b)  respec¬ 
tively.  For  the  coupling  coefficient,  we  used  the  following 
expression  from  Ref.  18, 


"SxV'V  . 


4/lM. 


(18a) 


V 


where  a 


e 


xy 


P 


M 

—  fl 

M,  L 


ce 

2 

°3 


(  1 


a)  a]  . 


(  18b) 


(  18c) 


We  have  defined  the  low  frequency  refractive  index  as 
n 3  =  k3c/w3  and  e 3X ,  e3y  are  the  elements  of  the  cold  plasma 
dielectric  tensor.20  Equations  (18a)  was  derived  under  the 
assumption  that  the  pump  and  its  sideband  are  unmagnetized 
(wi,2  >>  “’ce*  the  electron  cyclotron  frequency)  while 
the  low  frequency  wave  dispersion  relation  is  kept  qeneral 
but  with  w3  <  wpj.,  the  ion  plasma  frequency.  In  the  ELF 
range,  (u>3,  k3)  is  the  compressional  Alfven  branch  whereas  in 
the  VLF  regime,  it  corresponds  to  the  whistler  mode.  The 
calculations  were  performed  with  plasma  parameters  pertinent 
to  the  particular  altitude  under  daytime  conditions. 


The  results  for  VLF  generation  by  collisionless  inter¬ 
action  are  rather  disappointing  due  to  the  large  powers 
required  for  complete  decay.  Figure  4  shows  that  extremely 
high  electric  fields  are  needed  in  the  interaction  region  to 
overcome  the  decay  threshold.  If  no  resonant  enhancement  of 
the  field  occurs  in  the  interaction  region,  the  effective 
radiative  power  (ERP)  required  will  be  of  the  order  of  106- 
107  MW.  Such  a  peak  power  requirement  of  1-10  msec  duration 
is  beyond  the  capability  of  any  present  pulsed  power  systems. 
This  leaves  the  possibility  of  utilizing  the  near  field  of 
space-born  antenna  systems.  If  the  threshold  power  can  be 
achieved,  the  VLF  power  generated  will  be  of  the  order  of 
hundreds  of  KW. 


The  ELP  wave  generation  in  the  semi-coil isional  regime 
gives  more  interesting  results.  This  is  due  to  the  lower 
group  velocity  of  the  ELP  wave  and  the  collisional  reduction 
of  the  length  lc  as  given  in  (17c).  Figure  5  shows  a  much 
lower  threshold  field  which  lies  in  the  range  of  10-15  V/M 
for  a  frequency  of  around  100  Hz.  The  HP  ERP  in  this  case  is 
of  the  order  of  104  MW  and  for  times  less  than  1  sec,  this 
peak  power  requirement  is  well  within  the  state  of  the  art. 
The  threshold  can  be  further  reduced  if  the  increase  in 
ven  due  to  the  local  electron  HF  heating  is  included. 
The  downconverted  ELP  power  will  be  of  the  order  of  a  few 
hundred  watts  if  we  assume  an  antenna  gain  of  40  db.  It  is 
to  be  emphasized  that  the  above  values  are  for  the  specific 
case  of  HF  incidence  perpendicular  to  the  density  gradient  so 
that  a  key  requirement  for  a  proof  of  principle  experiment 
will  be  the  steerability  of  the  HF  antenna. 
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Three-wave-coupling  Interactions  have  been  examined  for  a  dissipative  medium.  It  is 
shown  that  in  the  semlcollisional  regime,  where  the  frequency  of  the  low-frequency  wave 
<„,)  is  smaller  than  the  collision  frequency  v  that  affects  the  weakly  damped  high-frequen¬ 
cy  waves  (a,,,  ^j),  i.e.,  >aj,  the  thermal  nonllnearitv  predominates  and  decay 

efficiencies  in  excess  of  that  expected  from  the  Manley-Rowe  relations  can  be  achieved. 

In  addition,  the  pump  threshold  power  for  complete  decay  is  substantially  reduced. 

PACS  numbers:  52.35.Mw 


The  interaction  of  two  high-frequency  wave 
packets,  U',, it,)  and  (w2,Ic2),  with  a  low-frequen¬ 
cy  one  (wj,ic3)  is  fundamental  in  the  study  of  non¬ 
linear  processes  in  dispersive  media.  An  often 
desired  goal  is  the  efficient  down-conversion  of 
high-frequency  power  to  low  frequencies.  Such  is 
the  case  for  rf  heating  of  plasmas,1  and  for  the 
development  of  wireless  ionospheric  communica¬ 
tion  schemes.2  For  interactions  in  nondissipative 
media,  the  maximum  energy  transfer  to  the  de¬ 
cay  products,  as  described  by  the  Manley-Rowe 
relations,3  is  proportional  to  their  frequencies, 
i'j'di,.  In  this  Letter  we  demonstrate  that  a  chan¬ 
nel  of  decay  exists  in  a  semicollisional  regime 
with  down-conversion  efficiency  in  excess  of  that 
predicted  by  the  Manley-Rowe  relations.  It  is  al¬ 
so  shown  that  there  is  a  significant  reduction  in 
the  threshold  required  for  complete  decay1  (i.e. , 
pump  depletion).  The  scaling  in  the  semicolli¬ 
sional  regime  is 

u)lfW2»»»U>3,  (1) 

where  v  is  the  collision  frequency  which  controls 
the  damping  of  the  high-frequency  waves.  In  this 
regime,  the  dominant  part  of  the  low-frequency 
nonlinear  force  on  the  electrons  is  due  to  the 
thermal  nonlinearity  (i.e. ,  the  same  nonlinearity 
that  produces  ionospheric  self -focusing*  or  en¬ 
tropy  modes6)  rather  than  the  usual  ponderomotive 
nonlinearity.  This  leads  to  nonsymmetrlc  coup¬ 
ling  coefficients  in  the  three-wave  interaction, 
with  very  profound  consequences.  Although  our 
results  are  general,  tor  concreteness  we  focus 
on  the  excitation  of  a  compressional  Alfv6n  wave 
by  a  high-frequency  pump  in  the  lower  ionosphere, 
where  the  required  scaling  (1)  is  satisfied.  This 
particular  Interaction  has  Important  applications 
to  wireless  communications  at  extremely  low 
frequencies. 

Consider  the  interaction  in  the  semlcollisional 
regime  described  by  (1),  with  v  given  by  electron- 


neutral  collisions.  For  the  high-frequency  pump 
and  its  sideband  we  assume  ,fl, , 

where  u>#  and  ft,  are  the  electron  plasma  and 
cyclotron  frequency.  The  low-frequency  wave  is 
taken  to  be  the  compressional  Alfv&n  branch. 
Following  standard  procedures,7  we  find  that  in 
a  homogeneous  medium,  the  slowly  varying  en¬ 
velopes  of  the  wave  packets  E ,  (S ,  t)  =  a1  a ,  (x ,  t ) 
*exp(-  t(w,<  - k; •  x)j,  subject  to  the  resonant 
conditions, 

k,  =  k j  +  k3,  (2) 


satisfy  the  system  of  equations 


3  „  3  _  ' 

.  V  \ 

3  _  3  _  ' 

K~*( 

,  V  \ 

3TV3i  +  r* 

'-fjvc. 

3  _  3  ' 

K-*( 

4  v  \ 

1+,3-3ha* 

(3) 


where  u,  are  the  group  velocities  and  T,  are  the 
damping  rates.  The  normalization  factors  are 
al  ,  =  c/VrA  and  o3  =  (w3/w,)l/2.  In  a  two-dimen¬ 
sional  geometry  (x  ,z)  with  the  ambient  magnetic 
field  B 0=Bgi,  K  is  given  by 


K=i 


JL  ka.  .  JfjL 
4m  (uijWj)1  2  u>,ft,  ’ 


(4) 


which  is  the  symmetric  coupling  coefficient  for 
the  collisionless  interaction  (u  =  0).  This  expres¬ 
sion  torK  is  valid  when  ^,<<0, ,  the  ion  cyclo¬ 
tron  frequency.  In  deriving  the  low-frequency 
equation  in  (3),  we  have  included  in  addition  to 
the  familiar  ponderomotive  nonlinearity,  the  non¬ 
linearity  due  to  Ohmic  heating  of  the  electrons,5 
which  produces  a  temperature  perturbation  T3  at 
u>j.  Using  the  electron  heating  equation*  we  ob¬ 
tain 


_  ,  4  v  e3E  r  £,* 

fa, - —  t  *  - 

3  3ul,  l»l 


(5) 
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where  the  electron  quiver  velocities  have  been 
approximated  by 


The  total  low-frequency  nonlinear  force  on  the 
electrons  Is  the  sum  of  the  thermal  force  due  to 
(5),  -n0VT3,  and  the  usual  ponderomotlve  force, 

f  pa-  V  fa*  +  V  $t) 

.£b(V|xbs*4-VxB1).  (?) 

where  n<,  is  the  ambient  plasma  density.  If  we  let 
v-*k  and  use  Faraday’s  law  for  t  In  (7),  the 
resulting  nonlinear  electric  field  Is  given  by 

■  +  — (8) 
5  JV  3uij  )m  Wji* 

which  then  drives  a  nonlinear  current  to  generate 
a,.  It  is  clear  from  (8)  that  the  thermal  nonlin¬ 
earity  dominates  when  v/u ).»  1,*  introducing  an 

...  ^  Mf 


wave  to  Wj/wj. 

Oir  analysis  shows  that  these  conclusions  are 
seriously  modified  when  one  considers  the  three- 
wave  interaction  in  the  semicollisional  regime. 
From  (3),  we  obtain  the  linear  equations  for  the 
eigenvalue  problem  that  determines  absolute  in¬ 
stability: 

(^+Ms,g7*r2  )oa  =  -y0(.t)o,*,  (11a) 

{h*u*h*T >)***“  (llb) 

where  y„  now  contains  the  pump  inhomogeneity 
and  A  is  defined  by 

A  *  1 +ff,v/,wJ#  (12) 

For  the  threshold  calculation,  the  small  colllsion- 
al  correction  on  the  right-hand  side  of  (11a)  has 
been  ignored.  If  we  use  Laplace  transformation 
in  time  and  make  the  substitution 


earuy  aommaies  wnen  >u  ,  .  ri/firnrN' 

equivalent  enhancement  factor  in  E  3nl  as  well  as  ^aa*,a3  j  =expj W  ***,), 

•>  „/2  nhosa  riifforonre.  Coimlinff  of  a.  to  the  hiah-  '  M**  Uix 


a  »/2  phase  difference.  Coupling  of  a,  to  the  high- 
frequency  waves  is  accomplished  by  the  beating 
of  with  the  low-frequency  electron  density 
perturbation 

(9) 

T,*  being  the  electron  contribution  to  the  cold- 
plasma  dielectric  tensor. 

In  the  collisloniess  limit  (v  -  0),  one  recovers 
from  (3)  the  standard  three-wave  Interactions 
which  have  been  investigated  by  the  inverse  scat¬ 
tering  transform  (1ST)  and  by  numerical  simula¬ 
tion.'*’  10  For  the  case  of  space-time  evolution  In 
one  dimension,  it  was  found  that  the  threshold 
for  pump  depletion  Is  equivalent  to  the  pump  val¬ 
ue  a3  required  to  achieve  an  absolute  instability 
within  a  width  L .  The  eigenvalue  problem  which 
determines  the  conditions  tor  absolute  instability 
In  a  rectangular  pump  is  well  known11  (l.e.,  the 
backward  wave  oscillator).  For  pump  inhomoge¬ 
neity  along  x,  the  instability  is  absolute  if 

unu3i  <  0»  (10a) 

L  >  ,  (10b) 

where  Le  •  Utau3l  1 1/3 /y0  is  the  critical  width, 
whiler0aKa1  is  the  uniform  medium  growth  rate. 
Nonlinear  saturation  of  the  absolute  Instability  by 
pump  depletion  as  suggested  by  soliton  solutions 
from  both  1ST  and  numerical  simulation  limits  the 
fractional  energy  transfer  to  the  low-frequency 


the  conjugates  of  (11a)  and  (lib)  become 


(£*’)■ 

(A-) 


where 


2 1  i<  2,  •’<  3l  r 


The  eigenvalue  or  growth  rate  p  may  in  general 
be  complex.  For  a  rectangular  pump,  (13)  read¬ 
ily  admits  solutions  in  terms  of  trigonometric 
functions  and  for  u2l  >  0,  «3l  <  0,  the  boundary 
conditions  are  that  a2  is  zero  at  -  L  1  while  a3 
is  zero  at  L  ,  2.  One  then  finds  the  dispersion  re¬ 
lation 

(e,e-nV)l,lcot(e,e  -hVV  *(^  (c)-n(c  =  o. 


where 

8=tan*Ht-v  w3),  (C  =  LC  |MlJ.  (16) 

We  have  solved  (15)  numerically  in  the  limit 
v  w3 »  1  so  that  8  =  « '2  and  |  A !  ~  v  w3.  The  real 
part  of  the  normalized  growth  rate  nle  for  the 
first  four  growing  modes  is  shown  in  Fig.  1  as  a 
function  of  the  normalized  width,  L  /!c .  The 
critical  width  threshold  now  becomes 

L  =>  1.81. .  (17) 
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FIG.  2.  Wave  profiles  at  t  =25,0  for  the  collislonal 
(solid  lines)  and  colllslooiess  (dotted  lines)  interactions. 


FIG.  1.  Real  part  of  nle  (nip  vs  L.'lc  <£  ' LJ  for  the 
colllslooal  (colllslonless)  Interaction.  The  dotted  lines 
refer  to  the  colllslonless  case. 

For  comparison,  the  results  for  the  collisionless 
case  (v  =  0;  9  =  0,  A  =  1 ,  lc  —Lc )  are  also  plotted. 
For  a  given  pump  width  L ,  the  value  of  al  re¬ 
quired  to  achieve  absolute  instability  in  the  semi- 
collislonal  regime  is  smaller  by  the  factor  le  Lc 
=  I A  |  'l/a.  This  implies  a  signif  icant  reduction  in 
pump-power  requirement  !at| J,  i.e. ,  of  the 
order  (v/w,)'*.  For  the  lower  ionosphere,  this 
is  more  than  two  orders  of  magnitude  reduction. 

Although  the  threshold  reduction  is  of  great 
practical  interest,  an  even  more  important  phys¬ 
ics  issue  arises  from  the  effect  of  nonsymmetric 
coupling  coefficients  on  the  down- conversion  ef¬ 
ficiency.  To  address  this,  the  fully  nonlinear 
evolution  of  (3)  is  required  so  that  pump  deple¬ 
tion  is  taken  into  account.  It  can  be  easily  shown 
that,  even  in  the  absence  of  damping  T,  =  0,  the 
Manley-Rowe  relations  no  longer  hold  for  this  in¬ 
teraction.12  The  sums  of  the  actions  are  not 
time  independent  and  the  total  energy  is  no  longer 
an  invariant  cf  the  system,  so  that 

3iSr/a<^f  i/A,  (18) 

where 

Sr-LS,,  (19) 

i- 1 

Consequently,  (3)  Is  not  amenable  to  analytical 
solutions  by  1ST  for  lack  of  conservation  laws. 
Instead,  we  resort  to  numerical  solutions.  For 
simplicity,  we  consider  the  model  problem  of  a 
one-dimensional  space-time  evolution  with  (*,.  =  0, 
'*2i3  ",4j« .  and  T,  =  0.  We  choose  wt  =  10.1,  w, 

=  10.0,  vs  1.2,  and  w,  =  0.1  in  accordance  with  (1). 
As  a  comparison,  the  collisionless  interaction 


(v  =  0)  is  also  examined.  The  initial  wave  packets 
are  rectangular  with  !a,|  /|oJ  =0.01,  i  u,|  =0.0. 
For  a  given  pump  width  L  ,  we  use  i  ati  =  1.0  for 
the  collisionless  case  and  |aJ=|A!  ‘1/J  =  0.25  for 
the  collislonal  case  so  that  the  respective  normal¬ 
ized  widths  L/Lc  and  L/l:  have  the  same  value; 
i.e.,  L  /Le  =L  /lc  =  3.2.  This  ensures  from  Fig.  1 
that  there  is  only  one  initial  growing  mode  in 
either  interaction. 

The  time-asymptotic  nonlinear  behavior,  at  f 
=  25.0,  is  shown  in  Fig.  2  where  pump  depletion 
is  seen  to  occur  in  both  interactions.  In  the  col¬ 
lisionless  case,  the  emergence  of  symmetric 
pulses  identifiable  as  solitons  is  as  predicted  by 
the  1ST  solutions.4  The  collisional  interaction, 
on  the  other  hand,  exhibits  nonsymmetric  pulses 
with  a3»a2  which  suggests  a  different  action- 
transfer  picture  than  that  described  by  the  Man¬ 
ley-Rowe  relations.  The  preferential  transfer  of 
action  to  a3  is  understandable  in  view  of  the  fac¬ 
tor  v  u.'3»  1  in  the  low-frequency  equation  as 
compared  to  the  small  corrections  r  w2  and  v 
in  the  equations  for  and  a  2,  respectively. 
Physically,  it  means  that  additional  w3  photons 
are  generated  at  a  large  rate  via  the  temperature 
perturbation  with  a  net  loss  of  energy  to  the  me¬ 
dium.  We  have  calculated  the  down-conversion 
efficiency  &3(t)  Sx(t  =  0)  for  this  interaction  and 
found  an  order-of- magnitude  enhancement  (~  20) 
over  the  collisionless  case.  The  energy  ratio 
S3(t)  i S,(f),  which  for  the  collisionless  interaction 
equals  the  frequency  ratio  *■,  w,  =  0.01,  has  a 
much  larger  value  (0.8)  for  the  collisional  case. 
The  above  considerations  are  valid  for  pump 
powers  such  that  the  nonlinear  interaction  time 
is  faster  than  the  electron  heating  time. 
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In  our  analysis,  we  have  neglected  the  role  o( 
parallel  heat  conduction  by  assuming  in  (5)  w3 
where  x  is  the  heat- conduction  coeffi¬ 
cient  Heat  conduction  can  be  included  by  simply 
replacing  w,  with  a>3*!fcj,2x  in  (3),  (5),  and  (8) 
and  the  definition  of  A  in  (12).  An  analytic  gener¬ 
alization  of  the  analysis  to  the  case  of  diffusive 
modes9  where  6„2x  >;>  ,  is  possible  by  rescaling 
the  equations  because  the  coupling  coefficients  in 
(3)  become  predominantly  real.  The  modified 
Manley-Rowe  relations  for  this  case  have  the  ap¬ 
proximate  form 


In  summary,  we  have  shown  that  much  lower 
pump  power  threshold  for  complete  decay  is  re¬ 
quired,  and  substantial  increase  in  down-conver¬ 
sion  efficiency  can  be  achieved  in  three-wave  in¬ 
teractions  in  the  semicollisional  regime  (v  »  u>3). 
This  result  can  have  very  profound  consequences 
for  schemes  requiring  efficient  high-frequency 
power  down-conversion,  such  as  wireless  iono¬ 
spheric  communciations,  low-frequency  rf  plasma 
heating,  etc.  Preliminary  estimates  indicate  that 
pump  power  of  a  few  megawatts  can  efficiently 
down-convert  high-frequency  power  to  extremely 
low-frequency  power  while  satisfying  the  theoret¬ 
ical  constraints.  These  will  be  discussed  in  fu¬ 
ture  publications. 
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